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ABSTRACT 

We present high resolution optical spectra obtained with the Ultraviolet and Visual Echelle 
Spectrograph (UVES) at the Very Large Telescope (VLT) and 21-cm absorption spectra 
obtained with the Giant Metrewave Radio Telescope (GMRT) and the Green Bank Tele- 
scope (GBT) of five quasars along the line of sight of which 21-cm absorption systems at 
1.17 < z < 1 .56 have been detected previously. We also present milliarcsec scale radio 
images of these quasars obtained with the Very Large Baseline Array (VLBA). We use the 
data on four of these systems to constrain the time variation of x = g p a 2 /fi where g p is the 
proton gyromagnetic factor, a is the fine structure constant, and yt is the proton-to-electron 
mass ratio. We carefully evaluate the systematic uncertainties in redshift measurements using 
cross-correlation analysis and repeated Voigt profile fitting. In two cases we also confirm our 
results by analysing optical spectra obtained with the Keck telescope. We find the weighted 
and the simple means of Ax/x to be respectively -(0.1 + 1.3) x 1(T 6 and (0.0 + 1.5) x 10~ 6 
at the mean redshift of < z >= 1.36 corresponding to a look back time of ~9 Gyr. This is 
the most stringent constraint ever obtained on Ax/x. If we only use the two systems towards 
quasars unresolved at milliarcsec scales, we get the simple mean of Ax/x = +(0.2+ 1.6)x 10~ 6 . 
Assuming constancy of other constants we get Aa/a = (0.0 + 0.8) x 10~ 6 which is a factor of 
two better than the best constraints obtained so far using the Many Multiplet Method. On the 
other hand assuming a and g p have not varied we derive Afi/fi = (0.0 + 1.5) x 10~ 6 which is 
again the best limit ever obtained on the variation of ft over this redshift range. Using inde- 
pendent constraints on Aa/a at z < 1.8 and Afi/fi at z ~ 0.7 available in the literature we get 
Ag p /g p < 3.5 x 10- 6 (lcr). 

Key words: galaxies: quasar: absorption line - quasar: individual: J0108-0037 - quasar: 
individual: J0501-0159 - quasar: individual: J1623+0718 - quasar: individual: J2340-0053 
- quasar: individual: J2358-1020 



1 INTRODUCTION 

Most of the successful physical theories rely on the constancy of 
a few fundamental quantities such as the fine structure constant, 
a = e 2 /he, the proton-to-electron mass ratio, fj, etc. Some modern 
theories of high energy physics that try to unify the fundamental 
interactions predict the variation of these di mensionless fundamen- 
tal constants over cosmological scales (see IUzanll2003L and refer- 
ences therein for more details). Current laboratory constraints ex- 



* Based on data obtained with UVES at the Very Large Telescope of the 
European Southern Observatory (Prgm. ID 082.A-0569A, 085.A-0258A, 
66.A-0624A, 68.A-0600A, 072.A-0346A, and 074.B-0358A) 



elude any significant variation of these constan ts over solar sys- 
tem scales and on geological time scales (see I Olive & Skillman 
l2004l : |Petrov et al.ll2006| ; |Rosenband et alj|2008h . It is not observa- 
tionally/experimentally excluded however that they could vary over 
cosmological scales. Therefore, constraining time and spatial vari- 
ations of fundamental constants of physics will have a great impact 
on understanding the true behavior of the nature. 



ISavedofil (l956) first pointed out the possibility of using red- 
shifted atomic lines from distant objects to test the evolution 
of dimensionless physical constants. Initial attempts in this field 
mainl y used the Alkali-Doublet (AD) method to constrain a varia- 
tion JSavedofl1ll95rj; iBahcall & Schmi"dlll96l Iwolfe et alj|l97d 
iLevshakovl 1 1994 fvarshalovich et al.l 1 19961: ICowie & Songaila 
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19951 : IVarshalovich et al.l l2000l : iMurphv et al]|200ll: IChand et al.l 



2005). While the AD method is simple, and least affected by sys- 



tematics related to ionization and chemical inhomogeneities, the 
limits achieved on Aa/a = (a z - a )/a are not usually stringent 
Q. The most precise value report ed to date using thi s method being 
Aa/a = -(0.02 ± 0.55) x 10" 5 dChand et alj|2005h over a redshift 
rang e of l.59<z< 2.92 

iDzuba et al.1 dl999dlbh and lWebb etai] d 19991) introduced the 
Many-Multiplet (MM) method as a generalization of the AD 
method, in which one correlates different multiplets from differ- 
ent ions simultaneously. Applying this method on a sample of 128 
absorb ers observed a t high spectral resolution with the Keck tele- 
scope, IMurphv et alj «{2003) claimed a detection, Aa/a=—(Q.51 ± 
0.10) x 10 , over the reds hift range 0.2< z <3 .7. However, this 
result was not confirmed by Srian and et al.l ( 12004b and lChand et al.l 
( 2004) who used higher signal-to-noise ratio (SNR ~ 70 per pixel), 
high spectral resolution (R >45000) UVES/VLT data of 23 Mgn 
systems detected towards 18 quasars in the redshift range 0.4 < 
z < 2.3 and fou nd Aa/a = (-0.06 ± 0.06) x 10~ 5 . This analysis 
was criticized by IMurphv et al.l d2007l) . However, from the reanal- 
ysis o f the UVES d a ta, us ing the Voigt profile fitting code VP- 
FItFI ISrianand et al. (2007) confirmed the null result albeit with 
larger error bars (i.e Aa/a = (0.01 ± 0.15) x 10~ 5 ). Other analysis 
using only Fen transitions in two particularly well suited absorp- 
tion systems at z = 1.15 and z = 1.84 failed to confirm any vari- 

i it 1 | - 1| — I | * 1 

ation in a (Ouast et al. 2004; Levshakov et al. 2006; Chand et al. 



l2006l ; lLevshakov et alj2007l) . Recently. I Webb et al.U201ll) have re 
ported the results of the analysis of 153 systems present in quasar 
spectra observed with VLT/UVES. They find that a increases 
with increasing cosmological distance fro m the Earth. Moreove r 
for z < 1.8, they confirm the results by Sria nand et al.l d2007l) . 
Aa/a = -(0.6 ± 1.6) x 10~ 6 . However combining their new VLT 
measurements with their previous Keck measurements they sug- 
gest the possibility for a spatial variation of a and speculate on 
the existence of an o-- dipole. If true, thi s dipole is very difficult 
to explain theoretically dOlive et ai1l201ll) . While the MM-method 
provides improved precision, it is affected by systematics related 
to ionization, chemical inhomogeneities and isotopic composition. 
The effects of inhomogeneities can be canceled using a large sam- 
ple of absorption systems but the effects of isotopic composition 
will likely to remain an issue. 

Most of the existing theories predict that the proton-to- 
electron m ass ratio shoul d vary much more than a (for ex- 
ample see lOlive et alj |2002| ; iDent & Fairbairnl 120031; | Dine et al.l 
I2003L and refere nces therein) though some predicts the reverse 
dPent et alj|2008h . The variations of fj. can be probed using H 2 Ly- 
man a nd Werner band absorption lines dVarshalovich & Levshakov! 
Il993h . H2 molecules are oc casionally dete cted in high redshift 
damped Lyman-q systems jpet itiean et al. 20 0j; ILedoux et aU 
120031 ; iNoterdaeme et al.||2008l ; ISrianand etal.ll2012l) with only a 
handful of them being suitable for probing the variation of /i. 
No clear indication of any variation in \i i n excess of 1 part in 
10 5 is seen in the ex i sting data for z > 2 dlvan chik et alj 2005 ; 
Reinhold et alj l200d: Ring et alj 120081; iThompson et alj 1 2009 ; 
Wendt & Molaroll201 ll ; lvan Weerdenburg et al.ll201 ll) . By compar- 
ing inversion line transitions of NH3 with the rotational transitions 
of other molecules, a strong constraint on Afi/fi can be obtained 



1 Here a z and ao are the measured values of a at any redshift, z, and in the 
laboratory on Earth 

2 http://www.ast.cam.ac.uk/ rfc/vpfit.html 



(Murphv et al. 2008). At present such an exercise is possible for 
only two g ravitationally lensed systems at z < 1 dHenkel et alj 
I2OO5U 2OO8). The b est reported c onstraint is Ayi/n < 3.6 x 10 7 (3cr) 
at z = 0.685 by iKanekan 1201 lh . Detecting more NH3 absorption 
towards normal quasars is require d to reduce systemati cs related to 
the usage of lensed quasars (See iHenkel et al. (2008) for discus- 
sions on various other systematics). 

As the energy of the 21 -cm transition is proportional to, 
x = a 2 g p //i, high resolution optical spectra and 21-cm spectra can 
be use d together to prob e the combined variation of these con- 
stants dWolfe et al.ll 19761) . Constraints of the o rder of cr(Ax/x) < 
10~ 5 were obtained towards individual systems j Cowie & S ongailal 
1 19951; iKanekar et alj|200l ; ISrianand et aT]|2010t) . iTzanavaris et alj 
d2007I) derived Ax/x = (0.63 ± 0.99) x 10~ 5 for a sample of nine 
21-cm absorbers with 0.23 < z < 2.35. The majority of the 
21-cm spectra used in this study were digitally scanned from the 
printed literature and the UV-optical data were obtained mainly 
with VLT/UVES. Better constraints can be derived from higher 
quality spectra in the radio and optical wavelength ranges of a well 
selected sample of 21-cm absorbers. This is possible now thanks 
to systematic surveys for 21-cm absorption towards strong Mg 11 
absorbers (e.g. lGupta e t al. 2009). This work has resulted in the de- 
tection of 9 new 21-cm absorption systems over a narrow redshift 
range (i.e 1.05< z < 1 .45) that can be used for constraining Ax/x. 

While this technique is very powerful there are two issues that 
introduce systematic uncertainties in the measurements. These are: 
(i) the identification of the optical component corresponding to the 
gas that produces the 21-cm absoiption and (ii) the fact that the 
radio and optical sources could probe different volumes of the ab- 
sorbing gas as the radio emitting region in quasars is in general 
extended compared to the UV emitting region. It has been sug- 
gested that the gas detected by their C i and/or H2 absorption is 
closely associated w ith the 21-cm gas dCowie & Songailal [19951 : 
ISrianand etall |2010|) . However, only few 21-cm absorbers show 
detectable C 1 and H2 absorption a nd even in these case s velocity 
offsets up to 1-2 km s -1 are noticed dSrianand et al.l2012l) . All these 
indicate C1/H2 and 21-cm absorption need not originate from the 
same physical region. Another option is to connect 21-cm absorp- 
tion to abso rption from singl y ionized species that trace Hi gas. 
For example lTzanavaris et al.1 fl2007h have associated the pixel with 
strongest absorption in the UV with the pixel with the strongest 
21-cm absorption. As neighboring pixels are correlated in optical 
spectra, the redshift of the strongest metal absorption component 
will be better defined by using simultaneous Voigt profile fits to 
the absoiption lines. This is the method we adopt in the analy- 
sis presented here. The second uncertainty discussed above can be 
minimized by selecting absorbers towards quasars that are compact 
at milliarcsec scales. While individual measurements may not be 
completely free of these systematics, even after careful considera- 
tion of the specific properties of the system, it should be possible 
to minimize them and get a statistically reliable measurement using 
large sample of absorbers. 

As different methods used for constraining the fundamental 
constants suffer from different systematic effects it is important to 
increase the number of measurements based on each method to ad- 
dress the time and space variation of different constants. Here we 
provide new measurements of Ax/x using a new sample of 21-cm 
absorbers. 

We have se lected 5 syste ms from the l iteratu re [4 from 
iGupta et al.1 d2009T) and one from lKanekar et al.1 d2009)] previously 
known to be associated with narrow 21-cm absorption lines towards 
radio sources that are compact at arcsecond scales. We have ob- 
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Table 1. Log of the optical spectroscopic observation with VLT/UVES 



Source name 


Exposure name 


Date 


Starting 


Exposure 


Setting 


Seeing 


Airmass 








Time (UT) 


(sec) 




(arcsec) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


J0108-0037 


EXP1 


2008-11-21 


01:15:02 


3700 


390+580 


0.64 


1.11 




EXP2 


2008-11-23 


02:41:45 


3700 


390+580 


0.77 


1.13 




EXP3 


2008-11-25 


01:52:08 


3700 


390+580 


0.89 


1.10 




EXP4 


2008-12-03 


01:11:20 


3690 


390+580 


0.71 


1.10 


J1623+0718 


EXP1 


2010-05-08 


05:48:49 


3340 


390+580 


1.33 


1.18 




EXP2 


2010-08-07 


01:08:41 


3340 


390+580 


0.82 


1.23 




EXP3 


2010-08-08 


00:45:37 


3340 


390+580 


0.86 


1.20 




EXP4 


2010-08-09 


00:56:50 


3340 


390+580 


0.70 


1.22 


J2340-0053 


EXP1 


2008-10-02 


02:39:30 


4500 


390+580 


0.88 


1.13 




EXP2 


2008-10-05 


02:01:02 


4500 


390+580 


0.82 


1.17 




EXP3 


2008-10-05 


03:25:10 


4500 


390+580 


0.88 


1.09 




EXP4 


2008-10-06 


00:28:59 


4500 


390+580 


0.90 


1.50 




EXP5 


2008-10-06 


01:55:32 


4500 


390+580 


0.79 


1.18 




EXP6 


2008-10-28 


04:41:03 


4500 


390+580 


0.79 


1.45 


J2358-1020 


EXP1 


2010-08-04 


06:27:07 


3340 


390+580 


0.75 


1.10 




EXP2 


2010-08-06 


04:43:45 


3340 


390+580 


0.69 


1.40 




EXP3 


2010-08-06 


05:50:01 


3340 


390+580 


0.71 


1.16 




EXP4 


2010-08-06 


06:55:07 


3340 


390+580 


0.64 


1.05 




EXP5 


7010-08-06 


0800- 14 

VJO . \J\J . L *T 


3340 


100+5R0 


0.65 


1.04 




EXP6 


2010-08-06 


09:05:21 


3340 


390+580 


0.77 


1.10 




EXP7 


2010-08-07 


04:53:34 


3340 


390+580 


0.78 


1.34 


J0501-0159 


EXP1 


2000-10-21 


06:13:17 


3600 


437+750 


0.61 


1.17 




EXP2 


2000-10-23 


04:08:46 


3600 


346+580 


0.53 


1.73 




EXP3 


2001-10-16 


07:20:51 


5400 


346+570 


0.46 


1.10 




EXP4 


2004-10-21 


04:38:08 


4500 


390+564 


0.63 


1.25 




EXP5 


2004-10-21 


05:42:38 


5400 


390+564 


0.79 


1.56 




EXP6 


2004-10-21 


07:05:51 


3600 


437+860 


1.17 


1.10 




EXP7 


2004-10-22 


04:38:08 


3600 


437+860 


0.84 


1.29 




EXP8 


2004-10-22 


05:42:42 


4500 


390+564 


0.89 


1.12 




EXP9 


2004-10-22 


07:05:55 


4500 


390+564 


0.59 


1.81 




EXP10 


2004-10-22 


08:04:58 


3360 


437+860 


1.00 


1.09 



Column 1: Source name; Column 2: Assigned name for the exposure. Column 3: Date of observation; Column 4: Starting time of exposure; Column 5: 
Exposure time; Column 6: Spectrograph settings; Column 7: Seeing in arcsec; Column 8: Airmass at the beginning of the exposures. 



tained high resolution UV and radio data of the quasars together 
with high resolution Very Large Baseline Array (VLBA) images. 
We report here the analysis of this dataset. This paper is organized 
as following. In Section 2, we present details of optical and radio 
observations and data reduction. In Sections 3 and 4 we provide de- 
tails of Gaussian fits to the 21 -cm absorption lines and Voigt profile 
fitting of the UV lines. In Section 5, we summarize our Ax/x mea- 
surements in individual systems and discuss the associated system- 
atic errors. In Section 6 we discuss the results and conclude. We 
use simultaneous Voigt profile fits to identify the redshift of the 
strongest UV component closest to the 21 -cm ab sorption. We also 
discus s the results if we adopt the method used bv lTzanavaris et ail 
d2007h . 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical spectroscopy 

The optical s pectroscopic observ ations of quasars were carried out 
with UVES flPekker et al1l2000h at the VLT UT2 8.2-m telescope 
at Paranal (Chile) in service mode [Programs 082.A-0569A and 
085.A-0258A], All observations were performed using the standard 
beam splitter with the dichroic #2 (setting 390+580) that covers 
roughly from 330 to 450 nm on the BLUE CCD and from 465 to 



578 nm and 583 to 680 nm on the two RED CCDs. Slit width of 
1 arcsec and CCD readout with 2x2 binning were used for all the 
observations resulting in a pixel size of ~ 1 .7 kms~' and spectral 
resolution of ^45000. 15"' Odorico et all d2000b have shown that the 
resetting of the grating between an object exposure and the ThAr 
calibration lamp exposure can result in an error of the order of a 
few hundred meters per second in the wavelength calibration. To 
minimize this effect each science exposure was followed imme- 
diately by an attached set of 5 ThAr lamp exposures. In the case 
of J0501-0159 we have retrieved all the UVES data available in 
the ESO archive. As these spectra were not acquired specifically 
for constraining the variation of fundamental constants there is no 
attached calibration lamp exposure taken along with the science 
observations. However this data were reduced using the available 
lamp spectra closest in time. 

The data were reduced with the UVES Common Pipeline Li- 
brary (CPL) data reduction pipeline release 4.7. fusing the optimal 
extraction method. We used 4 th order polynomials to find the dis- 
persion solution. The number of suitable ThAr lines used for wave- 
length calibration was always larger than 400 and the rms error was 
found to be in the range 70 - 80 m s -1 with zero average. However, 
this error applies only to regions very close to the ThAr emission 

3 http://www.eso.org/sci/facilities/paranal/instruments/uves/doc/ 
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lines that are used to compute the wavelength solution. In principle 
the calibration error in the regions in between ThAr emission line 
can be typically of the order of a few h undred meters per second 
(for example see lAgafonova et al]|201 ll) . 

All the spectra were corrected for the motion of the obser- 
vatory around the barycenter of the Sun-Earth system. The velocity 
component of the observatory's barycentric motion towards the line 
of sight to the object was calculated at the exposure mid point (see 
Table [TJ. Conversion of air t o vacu u m wa velengths was performed 
using the formula given in lEdlenl Jl966l) . For the co-addition of 
the different exposures, we interpolated the individual spectra and 
their errors to a common wavelength array and then computed the 
weighted mean using weights estimated from the errors in each 
pixel. In order to fit the continuum we considered only specific re- 
gions (20-100 A) around the absorption lines of interest and fitted 
the points without any absorption with a lower order cubic spline. 

Voigt profile fits of the absoiptions from different species have 
been performed using VPFIT, version 9.5. While simultaneously 
fitting absorption profiles of a system, we assumed that all the 
singly ionized species (e.g. Fen, Sin, Zmi, etc.) are kinematically 
associated with the same gas. We also assumed the velocity broad- 
ening is predominantly turbulent. Therefore, we used the same z 
and b parameters for a given component for all the species. The er- 
ror on the redshift of individual Voigt profile components depends 
on the statistical error from the fitting procedure and the systematic 
errors related to the procedure itself and the wavelength calibration. 
The VPFIT program estimates errors using only the diagonal terms 
of the covariance matrix. Although the reliability of the errors have 
been confirmed for unblended components (see iKing et alj|2009l ; 
ICarswell et al.l201 1), errors from VPFIT are underestimated in the 
case of blended components. To account for this and other system- 
atic errors discussed above, we perform Voigt profile fits for a given 
system several times (see Section |2.1.3| for details). In Table[2]we 
summarize the laboratory wavelengths and oscillator strengths of 
all transitions that are used in this study. In this table we also give 
a short name (id) to specific transitions for future reference in the 
text. 

2.1.1 Systematic errors in wavelength calibration 



The shortcomings of the ThAr wavelength calibration of quasars 
spectra taken with VLT/UVES has alr e ady been discussed by 



number of authors (Chandetal.1 20061 iLevshakov et al.l 


2006; 


Molaro et alj|2008l; iThomcson et al. 20091: Whitmore et all 


2010; 



Agafonova et al.l201 ll) . Overall velocity shifts of the order of a few 
hundred meters per second have been observed between the spectra 
of the same object calibr ated with an iodine ce ll spectrum or with 
a ThAr calibration lamp (Whitmore et al. 2010). The latter authors 
have also shown that intra-order velocity shifts of more than 200 
m s -1 are present within a given exposure. Therefore, while we have 
taken enough care during the data reduction, wavelength uncertain- 
ties due to these systematics still remain. We therefore performed 
several tests to estimate these systematic effects. 

2.1.2 Cross-correlation analysis 

We cross-correlate individual spectra (in a window comprising of 
the absorption lines associated with the systems of interest) with the 
combined spectrum to estimate the velocity offset between them. 
For this we first rebin each pixel (6v~2.0 kms~') to 25 sub-pixels 
with a velocity width of ~80 m s . This rebinning is done by inter- 
polating the spectrum with a polynomial using Neville algorithm 



Table 2. Adopted atomic data for different species used in this study 



Species id wavelength Ref oscillatorf 







(A) 




strength 


Sin 


a 


1808.01288 


[2] 


0.00208 


Cm 


bl 


2056.25682 


[2] 


0.1030 


Cm 


h2 


2062.23594 


[2] 


0.0759 


Cm 


b3 


2066.16391 


[2] 


0.0512 


Mn ii 


cl 


2576.87534 


[2] 


0.361 


Mn ii 


c2 


2594.49669 


[2] 


0.280 


Mini 


c3 


2606.45883 


[2] 


0.198 


Fe ii 


dl 


1608.45081 


[3] 


0.0577 


Fe ii 


d2 


1611.2005 


[1] 


0.00138 


Fe ii 


d3 


2249.8768 


[1] 


0.00182 


Fe ii 


d4 


2260.77934 


[21 


0.00244 


Fe ii 


d5 


2344.21282 


[21 


0.114 


Fe ii 


d6 


2374.46013 


[2] 


0.0313 


Fe ii 


d7 


2382.76411 


[2] 


0.320 


Fe ii 


d8 


2586.64937 


[2] 


0.0691 


Fe ii 


d9 


2600.17223 


[21 


0.239 


Nin 


el 


1454.842 


[1] 


0.0276 


Nin 


e2 


1467.259 


[1] 


0.0063 


Mn 


e3 


1467.756 


[1] 


0.0099 


Mn 


e4 


1502.148 


[1] 


0.006 


Mn 


e5 


1703.4119 


[1] 


0.006 


Mn 


e6 


1709.6042 


[1] 


0.0324 


Nin 


e7 


1741.5531 


[1] 


0.0427 


Mn 


e8 


1751.9157 


[1] 


0.0277 


Zmi 


fl 


2026.13695 


[2] 


0.501 


Zmi 


f2 


2062.66028 


[2] 


0.246 


Ci 


gl 


1560.3092 


[1] 


0.0774 


Ci 


g2 


1656.9284 


[1] 


0.149 


Mgi 


hi 


2026.47504 


[21 


0.113 


Mgi 


h2 


2852.96282 


[2] 


1.83 



References. Ill -iMortori 420031) ; [2] - lAldeniusj 12009]) ; [3] - 
iNave & Sansonettil 1201 ll) 

t Oscillator strengths are from Morton (2003) apart from the Niii/11454 
where we use the value from Zsargo & Federman ( 1998) that provide best 
fit to our data w hile being consistent within error to the values given in 
iMortorjfiool . 




4395.5 4396.0 4396.5 4397.0 4397.5 4398.0 
Wavelength (A) 



Figure 1. Example of the interpolation of a spectrum using Neville's algo- 
rithm. The histogram curve is an absorption profile and the smooth curve 
gives the interpolated spectrum achieved using Neville's algorithm. 
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combinations of spectra excluding one exposure at a time. This ex- 
ercise allows us to understand the influence of individual exposures 
on our final redshift measurement. Similarly using the final com- 
bined spectra we perform repeated Voigt profile fitting including 
and excluding different transitions. This will allow us to estimate 
the redshift uncertainties due to the choice of lines used in the Voigt 
profile fitting and also the random intra-order shifts. 

It is known that the wavelength calibration is most accurate 
in regions close to the ThAr emission lin es that are used to de - 
rive the pixel to wavelength solution (see Agaf onova et alj|201 lh . 
Therefore, we performed Voigt profile fitting of those transitions 
that have at least one ThAr emission line (that was used for wave- 
length calibration) within ±50 km s -1 from the UV-optical compo- 
nent that coincides with the 21 -cm component. 

We use the results of the above exercises along with the results 
of the cross-correlation analysis to quantify the final errors in the 
redshift measurements. 



Figure 2. Velocity shift between the combined spectra and one of the ex- 
posures for the Fen /11608 absorption profile in J2340-0057. The dashed 
curve shows the x 2 (see l2.1.lt as a function of the applied velocity shift and 
the continuous curve is the best fitted parabola. The short and long dashed 
vertical lines respectively mark the zero velocity and velocity of -315 m s~' 
at which the x 2 is minimum. 

(see Fig. QJ. For finding the relative shift between the combined 
spectrum / and the spectrum /j of the i'th exposure, we proceed as 
follows: we fix / and shift I, relative to / by steps of 0.04 times the 
original pixel size (i.e 80 ms~'). At each step characterized by the 
shift 8 A, a x 2 is calculated from the difference in fluxes between / 
and /, and the flux error of x 2 (SA) = Ej[7j - hj(SA)f /crUSA), 
where Zy is the normalized flux of the j'th pixel of the i'th expo- 
sure with the error o~u . This x 1 is a function of SA and is minimum 
when the two profiles are aligned. We fit the function x 2 (SA) with a 
parabola and the value of SA mm at which the^ 2 is minimum is taken 
as the wavelength offset between the two absorption profiles (See 
Fig. |2j. Following the standard statistical procedure we assign ltr 
errors to this <M ra j n by computing the required change in SA so that 
Ay 2 = x 2 — xL-.„ — 1 - This p rocedure is similar to that implemented 
by lAgafonova et al.l d20 1 lh to find the shift between their spectra 
but with the difference that they use the simple sum of the square 
of the differences in fluxes instead of^ 2 . Using x 1 has the advan- 
tage that we can associate an error to the measured shift. For each 
exposure we measure the shifts (of all transitions used in our Voigt 
profile fitting) along with the errors. Having measured a shift and 
error for each transition in each exposure we find their weighted 
mean and weighted standard deviation as an estimate of the sys- 
tematic error due to a constant shift in the redshift of the absorbing 
system. In addition, the correlation analysis not only allows us to 
identify exposures with abnormally large wavelength shifts but also 
to identify absorption lines that may be affected by calibration un- 
certainties in one of the exposures. Results of this exercise for four 
quasars in our sample are summarized in Tables IA1I IA2I I A3 1 and 
IA4I of the Appendix. It is clear from these tables that individual 
exposures have typical (rms) shifts of up to ~350 ms 4 . 



2.1.3 Repeated Voigt profile fitting analysis 

We use simultaneous fitting of several transitions to measure the 
redshift of a given component and the associated error. To estimate 
the latter we perform repeated Voigt profile fitting using several 



2.2 Archival Keck/HIRES spectrum 

For two quasars studied here (J2340-0053 and J0501-0159) high 
resolution echelle spectra were obtained with Keck/HIRES by 
Prof. Prochaska and collaborators a s part of their database archiv e 
for abundance studies in DLAs jProchaska et al.l 1200 ll |2007|) . 
The wavelength coverage in both cases is less than that of our 
VLT/UVES spectra. The spectral resolution of the J2340-0053 
Keck/HIRES spectrum is roughly the same as our VLT/UVES 
spectrum (i.e 6.0 km s ) but its SNR is less than ours. In the case of 
J0501-0159 the spectral resolution of Keck/HIRES is ~ 8.0 kms -1 
and both spectra have comparable SNR. We fit the absorption pro- 
files in the Keck/HIRES in order to compare the results obtained 
with the two telescopes. This exercise helps us to understand the 
systematic errors in the wavelength calibration and especially the 
existence of any global shift between the two spectra. 

2.3 GMRT observations and morphology of the background 
sources 

iGupta et all J2009T) used a bandwidth of 1 MHz split into 128 fre- 
quency channels in the course of their GMRT survey for 21 -cm 
absorption in strong Mg n absorbers. This yields a velocity reso- 
lution of ~4 kms~' per channel. In our new GMRT observations 
for three sources (J0108-0038, J2340-0053, and J2385-1020) we 
have used a band width of 0.25 MHz split into 128 channels yield- 
ing a channel resolution of ~ 1 km s~ 1 . To increase the spectral SNR, 
each object was observed for 24 to 30 hrs (i.e in three full synthesis 
observations). In the case of J1623+0718 a bandwidth of 0.5 MHz 
was used to adequately cover both the 21-cm components. Spectral 
resolution in this case is ~ 2 km s . 

The data were acquired in the two orthogonal polarization 
channels RR and LL. For the flux density/bandpass calibration of 
GMRT data, standard flux density calibrators were observed for 
10-15 min every two hours. A phase calibrator was also observed 
for 10 min every 45 min to get reliable phase solutions. The GMRT 
data were reduced using the NRAO AIPS package following the 
standard procedures. Special care was taken to exclude the base- 
lines and time stamps affected by the radio frequency interference 
(RFI). The spectra at the quasar positions were extracted from the 
RR and LL spectral cubes and compared for consistency. If nec- 
essary, a first-order cubic-spline was fitted to remove the residual 
continuum. The two polarization channels were then combined to 
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Table 3. Shifts (in units of m s 1 ) between individual 21 -cm spectra relative 
to the combined one for J0108-0037 at the position of 21-cm absorption 



LL[ RRi LLi 


RR 2 


LL 3 


RR 3 


47±55 99±96 71+104 


-180+133 


-178+191 


-249±192 


weighted mean 


+17 






weighted standard deviation 


122 







get the stokes I spectrum which was then shifted to the heliocentric 
frame. We used the AIPS task CVEL to correct the observed data 
for the Earth's motion and rotation. We obtained the mean spec- 
trum weighting the flux by the inverse rms square in the line free 
channels. 

2.3.1 Redshift uncertainties 

We can use the cross-correlation analysis described in Section |2.1.2| 
to search for any possible frequency off-set between spectra of the 
same object obtained at different epochs and through different po- 
larization channels. To avoid the effect of poor SNR, we consider 
here only J0108-0037. We have useful data obtained during three 
epochs and hence 6 individual spectra to carry out a correlation 
analysis of individual spectra relative to the final combined one. 
The results are summarized in Table [3] Here LLi and RRi are cor- 
responding to the LL and RR polarizations of the i'th observation 
of this source. The weighted standard deviation of these observed 
shifts is 122 m s -1 and observed values are up to 249 m s . We will 
thus consider the systematic error in radio frequency calibration to 
be 122 ms~' or equivalently 0.4xl(T 6 in Ax/x. 

Below, we will give the results of Gaussian fitting of the 21-cm 
lines using the combined spectrum and spectra obtained through in- 
dividual polarization channels. The statistical errors in 21-cm red- 
shift determination are found to be much larger than the above 
quoted systematic shift. 

2.3.2 Milliarcsec images of the background sources 

The five quasars being studied in this work seem to be compact in 
FIRST and GMRT 610-MHz images that both have 5 arcsec spatial 
resolution. VLBA L-band 1422 MHz observations of 4 sources (i.e. 
excluding J0501-0159) have been obtained as a part of a larger 
survey of radio sources with damped Lyman-ff absorbers (DLAs) 
and Mg n absorption systems along their line of sight to understand 
the relationship be tween radio structure and the detectabili ty of 21- 
cm absorption (see Srianand et al.l2012tlGupta et al .2012). Details 
of observational setups and data reduction can be found in these 
papers. The final images of these 4 sources are shown in Fig. [3] 

In Fig.[3]we see that J0108-0037 is clearly resolved into sev- 
eral components at milliarcsec scales. We note that 73% of the L- 
band flux density detected in the FIRST image is recovered in our 
VLBA image with the compact unresolved component having 53% 
of the flux density. In the 8 GHz VLBA image the strong compo- 
nent seen in our L-band VLBA image gets resolved into two dis- 
tinct components and only one of them may be associated to the 
optical continuum emitting region. This makes our task difficult 
while trying to assign 21-cm absorbers to their UV-optical absorb- 
ing counterparts as there is a high probability that additional con- 
tribution to the 21-cm absorption may come from gas that is not 
located along the optical line of sight. 

From Fig.[3] we see that the source J1623+0718 is unresolved 



even at milliarcsec scales. However, the VLBA observation recov- 
ers only 42% of the flux density detected in arcsec scale FIRST 
image. Based on a Gaussian fit to our VLBA image we derive that 
the maximum angular extent of the compact component is 4.64 
mas. This means the size of the radio beam at the redshift of the 
absorber, z ~ 1.337, is less than 39 pc. It is likely that this radio 
source samples the same region of the absorbing gas as the optical 
source. High frequency VLBA observations are not available in the 
literature for this source. In our recent L-band GMRT observations 
we find that the flux density of this source at arcsec scales is similar 
to what is seen in the FIRST image. This suggests that the ~58 % 
missing flux in our L-band VLBA image may be due to a diffuse 
extended radio emitting component that might have got resolved 
out at milliarcsec scale. Therefore if the absorbing gas is extended 
beyond 39 pc then the absorption against this diffuse component 
may also contribute to our GMRT spectrum. 

The source J2340-0053 is clearly unresolved in our L-band 
VLBA image. It is also unresolved in the h igh frequency VLBA 
images taken at 2 and 8 GHz (for example see lKovalev et alj2007t) . 
The radio flux density shows a peak around 1.4 to 2.3 MHz with a 
sharp decrease towards low frequency end. All this is consistent 
with the background quasar being a GHz peaked compact self- 
absorbed radio source. Using the flux density measurement in the 
FIRST catalog we find 90% of the L-band flux density in the FIRST 
image is recovered in the unresolved VLBA component. Gaussian 
fitting of the L-band VLBA image gives the maximum angular ex- 
tent of the object to be 1.71 mas. Using the redshift of the 21-cm 
absorber at z ~ 1.36 we estimate the maximum physical size of the 
quasar radio beam at the position of the absorbing gas to be < 15 
pc. Thus it is most likely that optical and radio beams sample the 
same volume of absorbing gas. Therefore we expect the system- 
atics related to the structure of the background radio source to be 
minimum in this case. 

The source J2358-1020, observed with a resolution of ~20 
mas (Fig.[3](, is unresolved in our L-band VLBA image. This source 
remains unresolved ev en in higher frequency VLBA observation s 
taken at 2 and 8 GHz dFomalont et"ai]|2OO0l : lFe"v & Charlotll2000l) . 
Similarly to the case of J2340-0053, this source is a GHz peaked 
radio source with a clear turnaround at the low frequency end. 
Moreover, in our VLBA observation we recover 74% of the flux 
detected in FIRST observation. A single Gaussian component fit 
represents its 1422 MHz VLBA image well. The largest angular 
size of the source is constrained to be <3.71 mas which is equal 
to <31 pc at the redshift of the absorbing system. All this suggests 
that the optical and radio sight lines probe the same volume of the 
absorbing gas. Therefore we expect the systematics related to the 
radio structure of this background quasar also be minimum. 

The background radio source J0501-0159, also known as 
PKS 0458-020, exhibits multiple components at arcsecond and 
milliarcsecond scale s. The radio s tructu re of this source is inves- 
tigated in detail by Brigg s et"ai] dl989h to determin e the spatial 
exten t of the z abs = 2.04 absorber 21-cm detected by (Wol fe et al.l 
Il985h . At 1.6 GHz, the source is resolved into two components 
that are mar ginally separated at 1 arcsec resolution (see Fig.l of 
iBriggs et al.ll 198911 . In the 10 mas resolution map at 608 MHz, the 
frequency that is also close to the redshifted 21-cm frequency of 
the z abs = 1.56 absorber, the compact 'core' at arcsecond scales is 
further resolved in to a jet and a diffuse component (see Fig. 2 of 
lBriggsetal.lll989h . The radio emission in these different compo- 
nents is strong enough to contribute to the detected 21-cm absorp- 
tion. Therefore, if the absorbing gas extends over several milliarc- 
secs, the possibility of velocity offsets between radio and optical 
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Figure 3. Contour plots of VLBA images at 1 .4 GHz. The restoring beam, shown as the ellipse, and the first contour level (CL) in mJy beam 1 are provided at 
the bottom of each image. The contour levels are plotted as CLx(-l, 1, 2, 4, 8,...) mJy beam -1 . 



absorption lines due to the extended radio structure cannot be ruled 
out in this case. 

3 GAUSSIAN FITS TO 21-CM ABSORPTION LINES 
3.1 J0108-0037 

Gupta et al. ( 2007) have reported 21-cm absorption from the z abs = 
1.371 Mgn system towards J0108-0037. Despite a low rest equiv- 
alent width of ~0.3 A for the Mgn doublet, other absorption lines 
of weaker metal transitions are clearly seen even in the SDSS spec- 
trum. Note that z a b s is very close to z cm with an apparent ejection 
velocity of only ~ 180 km s~ 1 . The 2 1 -cm absorption is well approx- 
imated by a single Gaussian component in the GMRT data with a 
resolution of 4 kms^'channeL 1 . Subsequently, we observed this 
system for 25 hours spread over three observing runs of similar du- 
rations to obtain a spectrum with 1 km s~' channel width. The final 
21-cm absorption spectrum is the inverse- variance weighted mean 
of these 3 spectra (see top panel in Fig. [4). It can be seen from the 
figure that the absoiption profile is smooth and that more than one 
Gaussian component is required to fit the profile. The smooth broad 
wing could indicate either a shallow component with a large veloc- 
ity width or a blend of several weak narrow components as sug- 
gested by the fit to metal absorption lines (see section l4~TV We find 



that at least three components are needed to obtain a good fit with a 
reduced x 2 of 1 . 1 . The result of this fit is over-plotted on top of the 
observed spectrum in Fig. [4] The redshift of the strongest 21-cm 
component is 1.3709710(11) (see Table[4]l. The reduced^ 2 is 1.03 
for a 4 component fit. This suggests that four components are ad- 
equate to represent the 21-cm profile. The redshift of the strongest 
21-cm component in this case is 1.3709694(53). 



3.2 J0501-0159 

The sight line towards this quasar is interesting as it covers two 21- 
cm absorbers at z = 2 04 and 1.56 dWolfe et alll985l : lKanekar et all 
2009). lKanekar et alj d2010l) fitted the 2 1 -cm absorption with a sin- 
gle component at z 2 i = 1.5605300(25) in a spectrum smoothed to a 
resolution of 1.3 km s^'channeL 1 . Using the pipeline based on the 
NRAO's GBTIDL package we re-reduced their arch ived GBT data. 
Details on the GBT data reduction can be found in lSrianand et al.l 
The original data was obtained with two sets of spectral 
resolutions (i.e. 0.33 and 0.66 kms~' channel -1 ). However, we re- 
binned the individual spectra to 1.3 kms~' channel -1 resolution 
before combining them. This is done to match the resolution to 
our GMRT spectra of other sources. The absorption profile of this 
21-cm absorber (shown in the top panel of Fig. |4) shows the ex- 
istence of two absorbing components. Redshifts obtained from the 



8 Rahmani, H. 



J01 08-0037, z Q = 1.370909 



0.00 



-0.05 




J0501-0159, z Q = 1.560535 



0.00 



-0.05 ■ 




-20 



20 



-20 



20 



J1 623+071 8, z Q = 1.335668 



0.00 



o 



6 -0.05 



-i — | — r-|— ■ — ■ — i — 1 — r — ■ — n — ' — ' — p — i — ' — r~ 
8 9 9 1 1 1 




=« t 

iJJ i Ljf 


cP ~X 

_1 1 ■ 1 ■ l_ 


H' 1/ ' ' ' ' 

llT i i i i 
iii i i 
iii i i 
iii i i 
iii i i 
i i. .i, . i. , . .i 



-20 



20 



40 



J2340-0053, z Q = 1.360857 



-0.5 



0.00 



-0.05 




-20 

Relative 



20 



velocity (km 







I 1 1 








i 1 1 
i 1 1 
i i j| 

1 h r 
i If 

i y 

i r 
. 1 1 








.173023 


20 



Figure 4. 21-cm absorption profiles of all the Mg 11 absorbers in our sample. The redshift used to define the zero velocity scale is given in each panel as zq. 
The ticks mark the positions of the 21-cm absorption components. The long-dashed vertical lines are the positions of UV-optical absorber components (see 
Section[4}. The normalized residuals (i.e. ([data] -[model]) / [error]) for each fit is shown in the top of each panel along with the lcr horizontal line. 



two components fit are summarized in Table [4] It can be seen that 
the spectra obtained in two polarization channels give consistent 
redshifts for the two components within measurement uncertain- 
ties (i.e 380 and 700 m s~' for the strong and the weak component 
respectively). 



3.3 J1623+0718 

The 21-cm absorption at z ~ 1 .336 towards J1623+0718 was first 
reported bv lGuptaet all li 2009) with two components. Using their 
spectrum we measure z 2 i = 1.3356755(53) and 1.3358518(108). 
As the background source as well as the absorption lines are weak 
we reobserved this system for 16.8 hours spread over two full syn- 
thesis with a channel width of 2 kms~'. Unfortunately the SNR 
in the final spectrum is not high enough to provide accurate red- 
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Table 4. Results of multi-component Gaussian fits to the 21 -cm absorption line 



Quasar 


Channel width 


Z21 


Sv 


SVRR 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 






1.3708647(115) 


1.45 






JO 108-0037 


1 


1.3709710(11) 
1.3710614(19) 


0.14 

0.24 






J0501-0159 


1 
1 


1.5605277(32)+ 
1.5605745(60) 


0.38 
0.70 


-0.12 ±0.16* 
+0.46 ± 0.67* 


+0.30 ±0.44** 
-0.16 ±0.40** 


J1623+0718 


4 


1.3356761(51)+ 


0.65 


+0.14 ± 1.42 


-0.47 ± 0.78 




1.3358591(98) 


1.26 


-1.60 + 3.70 


+0.76 + 2.41 


J2340-0053 


1 


1.3608595(14)+ 
1.3608874(106) 


0.18 
1.35 


+0.17 + 0.20 
+0.42 ± 0.48 


+0.00 + 0.19 
+0.39 ± 2.60 


J2358-1020 


1 


1.1730206(17) 
1.1730188(25)+ 


0.23 
0.35 


+0.74 + 0.51 
+0.61 ±0.48 


-0.55 + 0.28 
-0.30 ± 0.28 



Column 1 : Source name. Column 2: Channel width in km s . Column 3: Absorption redshift and its error. Column 4: The error in the absorption redshift in 

km s~' . Column 5: Measured velocity offset with respect to the zn given in Column 3 in kms~' when only RR spectrum considered. Column 6: Measured 

velocity offset in km s~' when only LL spectrum considered. 

* XX polarization in GBT; ** YY polarization in GBT 

+ component and corresponding Z2l used for measuring Ax/x 
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Figure 5. New and old combined 21-cm spectra of J1623+0718 respec- 
tively in panel (a) and (b) with the best fitted double Gaussian in each case 
shown as continuous curve. The thick histogram in the bottom of two panels 
show the residual of the fits that are shifted by a constant offset for clarity. 
The v = is at z = 1.335668. Although the redshift of the first component 
in both new and old spectra matches very well the redshift of the second 
component (at v ~ 35 kms~'in the top panel) differs by 7.7+2.0 kms~' 
between the two spectra. The normalized residuals (i.e. ([data]— [model]) / 
[error]) for each fit is shown in the top of each panel along with the lcr 
horizontal line. 



shallow feature. Because of this reason we do not use this compo- 
nent to constrain Ax/x. In order to increase the SNR further, we 
combi ned our new spectra with the spectra obtained by Gupta et al. 
(2009). The Gaussian fit to the combined spectrum is shown in Fig. 
[4] The fit results are summarized in Table |4] The measured z 2 j = 
1.3356761(51) for the main 21-cm component agrees well with the 
measurements based on spectra obtained in individual polarization 
channels. The redshift uncertainty in this case corresponds to a ve- 
locity of 650 ms~'. 

3.4 J2340-0053 

The 21-cm absorption at z = 1.360 6 along the line o f sight to- 
wards this quasar was discovered by Gupta et al. (2009). We have 
acquired additional 3x8 hours of GMRT data at 1 km s~' channel -1 
resolution. However these data are found to be unusable due to RFI. 
So we use here only the the GMRT 21-c m absorption spe ctrum 
with, Sv ~ 1 km s~' channel -1 , obtained by Gupt a et all d2009h and 
shown in Fig. [4] The spectrum clearly shows two absorbing com- 
ponents, the bluer being quite strong. The continuous line in Fig|4] 
shows the double Gaussian fit to the 21-cm absorption profile. The 
two components are also shown (see also Table[4}. The redshifts of 
the two components are z= 1.3608595(14) and z= 1.3608874(106). 
The redshift errors for the Gaussian fits correspond to uncertainties 
in the velocity scale of 180 ms~' and 1350 ms~' respectively. The 
above redshifts are found to be consistent with those obtained by 
fitting RR and LL spectra separately (see TableQ. 



shift measurements. Therefore, we rebinned our new spectrum to a 
resolution of 4 kms~' channel -1 . The redshifts of the two compo- 
nents in the new combined spectrum are Z21 = 1.3356782(85) and 
1.3359118(116). While the redshift of the blue compo nent is con- 
sistent with the measurement based on the spectrum of Gupta et al. 
(2009) the redshift of the second component is off by 7.7±2.0 
kms -1 (see Fig. [5}. We attribute this to the low optical depth in 
this component or to the presence of a low level RFI affecting the 



3.5 J2358-1020 

The 21-cm absoiption along the line of s ight towards this quasar 
was first discovered bv lGupta et al.l 120071) . The 21-cm absorption 
is well approximated by a single Gaussian compon ent at a spectral 
resolution of 2 km s -1 channel 1 jGupta et al.ll2009h . Two high res- 
olution spectra of this object with Sv ~ 1 kms -1 per channel were 
acquired during subsequent observations (2x8 hours). The shape 
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Figure 6. Voigt profile fits to the absorption profiles of the z,nbs ~ 1 37 system towards J0108— 0037. The histogram plot in each panel shows the observed 
absorption profile of a given transition and the continuous curve is the best Voigt profile fit (or Gaussian fit in the case of 21 -cm profiles). Last panel shows the 
21-cm absorption profile along with the fitted model and its individual components. The normalized residuals (i.e. ([data] -[model]) / [flux error]) for each fit 
is shown in the top of each panel along with the icr horizontal line. 

Upper vertical tick marks indicate the position of 21-cm absorber components and lower vertical tick marks indicate different optical-UV velocity 

component. 



of the 21-cm absoiption feature even at this higher resolution is 
consistent with a single component (see Fig. [4j» . This is the sim- 
plest profile in our sample. It is best fitted with a single component 
at zji = 1.1730206(17). The typical error in the redshift measure- 
ment is ~230 m s . As there are 4 individual spectra we also per- 
formed a fit using errors that are the rms of the fluxes measured in 
the different spectra. The fit obtained using these errors gives z%\ = 
1.1730212(20) which is consistent with the above quoted measure- 
ment. In Table [4] we also provide results of independent fits to RR 
and LL spectra. The zn measurement based on the RR spectrum is 
higher than the one obtained from the LL spectrum with a relative 
off-set of 1.3±0.6 kms~'. We find this is mainly due to one of the 
RR spectra being affected by low level RFI. While this spectrum 
does not influence the weighted mean I-spectra, the combined RR 
spectrum is appreciably affected by this. Combining all the spectra 



but this affected spectrum yields z%\ = 1.1730188(25) (i.e with a 
redshift error corresponding to 344 m s -1 ). This is the z 2 i we use to 
derive Ax/x for this system. 



4 VOIGT PROFILE FITTING OF UV LINES 

In this section we describe the Voigt profile fitting of the UV ab- 
sorption lines and discuss the individual systems in detail. 



4.1 System at z abs ~ 1.37 towards J0108-0037 

J0108-0037 is one of the brightest quasar in our sample with SDSS 
r-band magnitude of 17.5. The SNR in the continuum close to the 
absorption lines used for measuring redshifts is usually larger than 
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Figure 7. Same as Fig.|6]for z-jbs ~ 133 Mgii system towards J1623+0718. 



30 for this system. The absoiption profiles of all singly ionized 
transitions used in the Voigt profile fits are shown in Fig. [6] In- 
terestingly many transitions of Ni n are detected in this system. 
This could allow one to measure Aa/a by using only Nin tran- 
sitions once accurate values of rest frame wavelengths, oscilla- 
tor strengths, and sensitivity coefficients are available. Apart from 
Niii/11467, 1502, 1703, other absorption profiles are very strong. 
Absorption profiles of Fe nA 1608,2344,2374,2382,2586,2600 in 
our final combined UVES spectra are highly saturated and have not 
been used for redshift measurement. We have fitted this system with 
two, three, four, and five components to find the optimal fit. The re- 
duced^ 2 are respectively, 1.60, 1.42, 1.28, and 1.23. Increasing the 
number of components does not lead to any better fit. Therefore we 
consider the fit with five components as the best fit for this absorb- 
ing system. From our best fit, there are two strong UV absorption 
components (at v = 6.4 and 19.5 kms~' in Fig. [6j with approxi- 
mately the same column density of metals. This means, unlike in 
other cases discussed here, a unique identification of the stongest 
UV component to be associated with the strongest 21 cm compo- 



nent is highly questionable. Even though both 21-cm and UV ab- 
sorption lines span the same velocity range there is no one to one 
correspondence between the two. This could mean that the 21-cm 
optical depth does not scale with the column density of metal lines. 
To illustrate this we plot in the bottom right panel of Fig. [6] a five 
component fit of the 21 cm absorption profile. This could mean as 
well that additional contribution to 21-cm absorption comes from 
gas that is not probed by the optical sight line. As the morphology 
of the background source is complex, we can not rule out that the 
differences in the absorption profiles are due to the fact that the op- 
tical and radio sight lines probe different volumes of the absorbing 
gas. Therefore, because of the degeneracy introduced by this pecu- 
liar profile and the complex morphology of the radio emission we 
do not use this system for Ax/x measurements. 



4.2 System at z ~ 1.33 towards J1623+0718 

J1623+0718 with SDSS r-band magnitude of -17.5 is another 
bright quasar in our sample. The velocity plots of some of the 
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Figure 8. Same as Fig.[6]for the z a bs~l-36 Mgn absorber towards J2340-0053 and for the VLT/UVES. In the bottom right comer panel we zoom over 40 
km s to show the 21 -cm absorber with the two components fit overplotted. 



species detected from the z a t, s ~1.33 system are shown in Fig.[7] The 
SNR in the continuum close to absorption profile of Niii/!1454 is 
~10, and it is higher than 30 close to Feii/12586. The absorption 
profiles of weak metal transitions like Siii/11808 are spread over 
80 kms~'. The 21 -cm absorption including the broad component 
is spread over the same velocity range. From visual inspection it is 
clear that the strongest metal absorption component coincides well 
with the main narrow 21 -cm component. 

The optical absorption profiles suggest the presence of addi- 
tional weak components in the wings (at v~53 and — 20 kms^ 1 ). 
We constrain the component structure in the wings from the strong 
Fe n lines for which the wings are apparent. The best fitted Voigt 
profile shown in Fig.|7]has a reduced^ 2 of 1.04. The measured red- 
shift of the strongest metal component in the 21 -cm velocity range 



is Zabs = 1.3356684(19) with the typical redshift error of 240 m s -1 
(see Table [5} . 



As can be seen, the main component in the metal absorption 
is broad and may contain additional hidden narrow components. 
Therefore, we repeated the fits with additional components injected 
around the main component. The reduced x 1 does not change with 
the addition of these new components. There is however a minor 
change in the absorption redshifts, albeit with increased errors, for 
the strongest component we find z = 1.3356675(47) with an error 
of ~600 ms~'. This is the UV absorption redshift and associated 
uncertainty we consider for Ax/x measurement. 
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Figure 9. Same as Fig.|6]for the z abs ~1.36 Mgn absorber towards J2340-0053 and for Keck/HIRES. 



4.3 System at z abs ~1.36 towards J2340-0053 

In addition to our VLT/UVES spectrum we have also analysed the 
Keck/HIRES spectrum of this object. The velocity plot of different 
ions detected in this system is shown in Fig.[8]and Fig.[9] The SNR 
in the UVES continuum is ~40 close to Ni H/11709 and can be larger 
than 60 close to Mn n absorptions. The absorption lines of singly 
ionized species are spread over ~150 km s . Unlike in the case of 
J1623+0718 and J0108-0037 the 21 -cm absoiption is very much 
na rrower than the UV absorption lines. However, as pointed out 
by I Gupta et"al] {2009), the metal component associated with the 
21 -cm component is well detached. In particular the component is 
well defined by Siii/11808 and Znn lines. Note, we have not fitted 
the absorption profile of Sin for v < -55 kms -1 as this region is 
contaminated by absorption from another intervening system. The 



fit has a reduced^ 2 of 1.3. The overall profile is fitted with 14 Voigt 
profile components and the absorption coinciding with the 21 -cm 
absorption requires two narrow components. 



The redshifts of the two UV-optical components in our UVES 
spectrum that are closer to 21-cm absorption are 1.3608565( 16) and 
1.3608781(19). The former happens to be the stronger absorption 
component to be associated to the stronger 21-cm absorption com- 
ponent. The redshift errors due to Voigt profile fitting correspond 
to 200 and 240 m s for these components. The redshifts measured 
from the HIRES spectrum, 1.3608572(36) and 1.3608759(40), are 
consistent with those derived from the UVES spectrum. The fitting 
errors in the redshifts for the HIRES data are 457 and 508 ms~' 
respectively. 
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Figure 10. Same as Fig.|6]for the z a bs~l-17 Mgn absorber towards J2358-1020. 



4.4 System at Zlibs ~ 1.17 towards J2358-1020 

J2358-1020 with an r-band SDSS magnitude of 18.7 is one of the 
faintest quasars in our sample. The SNR in the UVES spectrum 
is at best <25. As it can be seen in Fig. [10] the absorption pro- 
files are spanning more than 150 kms -1 . Absorption profiles of 
Crn,l2056,2062,2066, Znnvl2062, Mn 1^2594,2606 are not used 
in the fit as they are weak or located in regions of poor SNR. As 
there are weak components at v — 140 and 15 km s _1 , we included 
two core-saturated profiles, Fen/12344 and Feii/12586, to be able 
to fit the overall profile. Similar to the case of J2340-0053 the 
metal component coinciding with the 21 -cm absorption seems to 
be well detached from other components and is well fitted with 
a single Voigt profile component. This is apparent for Siii/11808, 
Znn/12026, and Mn n/12576. The overall fit has 11 Voigt profile 
components with a reduced^ 2 of 1.2. 

The redshift of the stronger absorption component that we as- 
sociate with the 21-cm absorption is 1.1730227(29). The redshift 
error from line fitting is ~400 ms~'. 



4.5 System at z abs ~ 1.56 towards J0501-0159 

J0501-0159 is a faint quasar with an r-band magnitude of 19.33. 
Similar to the case of J2340-0053 we have spectra obtained with 
VLT/UVES as well as Keck/HIRES. Fig. [TJJ presents the velocity 
plot of different absorption profiles detected in this system. The top 
and bottom panels show the VLT/UVES and Keck/HIRES spectra 
respectively, along with their best fit profiles. Although the final 
combined UVES/VLT spectrum is made of 10 exposures each of 
more than 3300 s of exposure time, the typical SNR is only <20. 
Apart from those Fe n lines that are shown in Fig. [TT] other Fe n 
lines are highly saturated or have poor SNR. Associated Ni n ab- 



sorption lines are very weak and are not used for the Voigt profile 
fitting. Some of them are located in the Lyman-a forest, others are 
either contaminated or have poor SNR and therefore have not been 
included in the fit. Our best fit is shown in Fig.QjTand has a reduced 
X 1 of 1.0. 

It can be seen in Fig. [TJJ that the strongest 21-cm component 
coincides with the strongest metal component seen in the absorp- 
tion profiles of undepleted species like Zn n and Si n. The redshifts 
of the UV-optical component in UVES and HIRES data are, respec- 
tively, 1.5605354(43) and 1.5605398(276) with errors of 503 ms 1 
and 3234 ms~'. The two measurements agree within l.Ocr. 



5 CONSTRAINING Ax/x 

In this section we present Ax/x measurements for individual sys- 
tems and discuss the associated errors in detail. We measure Ax/x 
using, 



Ax 

x 



ZUV ~~ Z21 
1 +Z21 



(1) 



where zuv is the redshift of the UV-optical component as given in 
column 2 of Tableland Zii is the redshift of the the 21-cm compo- 
nent as given in column 3 of Table [4] Usually the metal absorption 
has a larger velocity spread compared to that of 21-cm absorption. 
So to associate the UV-optical absorption (i.e. zv\) to the 21-cm 
component we use the stronger absorption component in the ve- 
locity range of 21-cm absorption. In all cases discussed here the 
UV-optical absorption clump associated to the 21-cm absorption 
is easily identifiable. The statistical error due to the fits, c(— ), is 
calculated as: 
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Figure 11. Same as Fig.[5]for the z abs ~1.56 Mgn absorber towards J0501-0159 from VLT/UVES spectrum (top) and HIRES/Keck spectrum (bottom). 
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Table 5. Results of repeated Voigt profile fitting analysis. 



Quasar 

(i) 


zuvt 
(2) 


Sv 
(kms -1 ) 
(3) 


6 V \ 
(kms -1 ) 
(4) 


id 
(5) 


si 

(kms -1 ) 
(6) 


id 

(7) 


H 

(kms -1 ) 
(8) 


id 

(9) 


*! 

(kms -1 ) 
(10) 


Sv] 
(kms -1 ) 
(11) 


J1623+0718 


1.3356684(19) 


0.24 


-0.02+0.23 


d2,el,e5 


+0.34+0.31 


d5 


+0.36+0.26 


d2,d3,d7,d8,e5,c2 




+0.23+0.22 


J2340-0053 


1.3608565(16) 


0.20 


+0.26±0.18 


d2,fl 


-0.22+0.19 


d8 


-0.48+0.33 


e7,a,d3,cl 


+0.09+0.46 


-0.07+0.15 


J2358-1020 


1.1730227(29) 


0.40 


+0.26±0.47 


fl,d4 


-0.01+0.37 


d9 


-0.54+0.44 


d5,d6 




-0.02+0.23 


J0501-0159 


1.5605354(43) 


0.50 


+0.35+0.27 


b2,b3,d4,c3 


+0.26+0.30 


d9 


+0.15+0.67 


b3,a,d3,d8,c2 


-0.52+3.23 


+0.07+0.17 



Column 1: Source name. Column 2: The absorption redshifts and associated error (in brackets) measured using our VLT/UVES spectrum. Column 3: error in 
the redshift measurement given in Column 2 in kms -1 . Column 4: measured velocity offset when excluding weak transitions (as defined in Table.[2j listed in 
Column 5. Column 6: measured velocity offset when a saturated line (with ids given in column 7) is included in the fit. Column 8: Velocity offset measured 
after excluding the absorption lines (whose ids are given in column 9) far away from the ThAr lamp lines used for wavelength calibration. Column 10: 
Measured velocity offset for the redshift measured using Keck/HIRES spectrum. Column 1 1 : Mean and standard deviation of measured redshifts after 
removing one exposure from combined spectra (see Fig. 1121 . 

' all the velocities and associated errors are calculated with respect to the main redshift given in the second column. 



J0501-0159 1.5605354 



500 


-500 



1 2 3 

J1 623+071 8 1.3356684 



o 
_o 

CD 
> 



CD 



1 2 3 

J2340-0053 1 .3608565 



1 2 3 4 5 6 

J2358-1020 1.1730227 



500 


-500 



l:±: F t: 



2 3 4 5 6 7 

Excluded exposures 



Figure 12. Results of an experiment consisting of measuring the shifts in 
the observed redshift as a consequence of retrieving one exposure when 
combining the different exposures to obtain the final quasar spectrum. We 
only show those components that are used for measuring Ax/x. Abscissa 
in each panel indicates the missing exposure(s). The zero of the y-axis is 
taken at the redshift (indicated at the top of each panel) measured in the 
spectrum obtained by combining all exposures (dotted line). Short-dashed 
lines shows the lcr error on this redshift. 



cr( — ) = 

X 1 + Z21 



1 + Zy\ 
1 +Z21 



X erf + erf 



(2) 



in which cr_ uv and cr Z2l are the errors on z uv and z 2 i respectively. In 
the case of <r, 21 we consider the contributions of statistical (Table. 
[4]l and systematic uncertainties (Table |3}, the latter being smaller 
(~ 122 ms 4 ) than the former (> 300 ms 4 ), To estimate the sys- 
tematic errors in the optical redshifts we carry out a number of tests 
whose results are summarized in Table [5] and also in the tables of 
the Appendix (see Section [2~T1 for more detail). In columns (4), (6), 
and (8) of Table [5] we give velocity offsets measured with respect 
to z abs given in column (2) using repeated Voigt profile fitting after, 
respectively, excluding weak lines, including saturated lines, and 
excluding the absorption lines with no ThAr line within 50 km s _I 
of that optical component assigned to 21 -cm component. The refer- 
ence to the corresponding lines are given in the preceding column. 
The results of these tests are sensitive to the infra-order wavelength 
calibration errors. To be on the conservative side we consider the 
maximum error found here as a measure of the systematic error 
introduced from intra-order shifts (cr c ). 

In column (11) of Table [5] we give the mean velocity offset 
found by repeated Voigt profile fitting of lines after excluding one 
of the exposures (See also Fig. [T2t . Tables in the Appendix also 
summarize the results of cross-correlation analysis between the in- 
dividual exposures and the combined spectra. The results of the last 
two exercises (exposure removal and cross-correlation) are sensi- 
tive to any constant shift between different exposures. Therefore 
we consider the maximum of these two shifts as an estimate of the 
constant shift error in the wavelength calibration (cr c ). The final 
systematic error for each sight line is taken as the quadratic sum of 
the two systematic errors. The approach taken here is conservative 
to allow for maximum uncertainty in the measurement of zuv- We 
summarize the values of <x e and cr c for individual systems in Table 
[6]where these errors are converted from velocity shifts to Ax/x. In 
column (7) of Table [6] we present the total systematic error, <r sys , 
which is calculated from the quadratic sum of cr e , cr c and 122 m S" 1 
we found from 21 -cm analysis. 



6 RESULTS AND CONCLUSIONS 



In Table[6]we summarize the Ax/x measurements in individual sys- 
tems. We recollect that the values are obtained under the assump- 
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Table 6. Ax/x measured (in units of 10 6 ) from different absorption systems 



Quasar 


^■abs 








UVES 






Keck 




dipole 


A(Aa/a) 






Ax/x 


0"stal 


<r c 


<r e 


0"sys 


0"tot 


Aa/a 


Ax/x 





Aa/a 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(ID 


(12) 


(13) 


J1623+0718 


1.3356 


-3.7 


3.0 


1.0 


1.2 


1.6 


3.4 


-1.8±1.7 




66.0 


+3.9+1.6 


+5.7+2.3 


J2340-0053 


1.3608 


-1.3 


0.9 


0.7 


1.6 


1.7 


2.0 


-0.6±1.0 


-1.0±1.6 


90.6 


-0.1 + 1.2 


+0.5+1.6 


J2358-1020 


1.1730 


+ 1.8 


1.8 


0.8 


1.8 


2.0 


2.7 


+0.9±1.4 




84.9 


+0.8+1.1 


-0.1 + 1.8 


J0501-0159 


1.5605 


+3.0 


2.1 


0.6 


2.2 


2.3 


3.1 


+ 1.5±1.6 


+4.7±10.9 


119.0 


-5.2+1.8 


-6.7+2.4 



UVES Keck 
simple average 0.0+1.5 +1.8+2.8 

weighted average -0.1+1.3 -0.9+1.6 



Column 1: Object name. Column 2: absorption redshift. Columns 3-8: Ax/x measurement, associated statistical and systematic errors (as discussed in Section 
5) respectively. Column 9: Aa/a calculated based on the final values for Ax/x assuming constancy of other constants. Column 10: Ax/x measu rements based 
on Keck/HIRES data. Column 11 and 12: angular distance in degrees between the quasar sight line and the best fitted dipole position from Web b et all ilOl ll) 
and the predicted value for Aa/a based on dipole. Column 13: difference between our measurement and the prediction from dipole and its associated error. 
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Figure 13. Comparison of our Ax/x measurements with three other mea- 
surements in th e literature. Filled circ les are our measurements and the filled 
square is from Srianand et al. (2010) using a 21-cm absorber at z ~ 3.174 
along the line of sight of J1337+3152. The dashed-dotted box corresponds 
to the standard deviation of our four measurements The l<r error around 
the mean is shown as a solid box. T he dashed line and long -dashed box are 
the Ax/x and its error measured by Tzanavaris et al lj2007l) . The filled and 
empty stars are respectively from iKanekareTalJ l l2010l) "and lKanekar et alj 
(2006). 



tion that the strongest UV and 21-cm absorption are produced by 
the same gas. The final errors in Ax/x for our VLT/UVES mea- 
surements given in column (8) are the quadratic sum of the statis- 
tical and systematic errors. We find the simple mean of Ax/x re- 
gardless of the associated errors in individual measurements to be 
(0.0±1.5)xlfT 6 with an rms of 3.0xl0~ 6 around the mean. A con- 
stant Ax/x of O.OxlfT 6 has a reduced^ 2 of 1.0 for our four UVES 
measurements that shows the estimated errors in Ax/x are not un- 
derestimated. If we apply the standard procedure of weighting the 
data points by their inverse square errors (given in column (8) of 
Table|6]l, we get -(0.1 ± 1.3) x lO" 6 . 

Our VLBA images suggest that two of the quasars in the sam- 
ple (i.e J1623+0718 and J0501-0159) may have resolved struc- 
tures at milliarcsec scale containing more than 50% of the flux. 



This may imply that if the absorbing gas is extended then some 
additional 21-cm absorption can originate from the gas that is not 
probed by the optical sight lines. However in the remaining two 
quasars this is not the case as most of the flux is recovered in the 
unresolved VLBA component. If we only use these two cases we 
derive Ax/x = +(0.2± 1.6)x 10~ 6 . This is very much consistent with 
what we find using all the four systems. Therefore the analysis pre- 
sented here does not find any statistically significant variation in x 
and this null result may not be related to systematics due to radio 
structure. 

There are two quasars for which we have spectra from both 
VLT and Keck. As can be seen from Table [6] in both cases the 
VLT and Keck measurements are consistent with each other with 
Keck measurements having larger statistical uncertainties. The 
mean Ax/x from VLT/UVES data of (0.0±1.5)xl0~ 6 is consistent 
with Ax/x = +(1.8±2.8)xl0~ 6 from Keck/HIRES data where we 
could not include the systematic error. This is also the case for the 
weighted means. 

In Fig.[T3]we compare our Ax/x results with other Ax/x mea- 
surements from the literature. Filled circles are our measurements. 
The dashed-dotted rectangular box indicates the mean and standard 
deviation of our measurements. The solid rectangular box gives the 
mean and final error on it when combining the 4 measurements. 
The green b ox with the dashed line gives the weighted mean and 
lcr found bvlTzanavaris et all d2007l) . The two filled stars are from 
iKanekar et alj l201Cl) and the empty star is from iKanekar et al.l 
(2006). The data point at z ~ 3.174 is from lSrianand et al.l d201oll 

The better accuracy reached in our study is mainly due to the 
following reasons: (1) Systems are chosen to have narrow 21-cm 
absorption components. (2) Three of the quasars have high resolu- 
tion (R~45000) and high SNR UV-optical spectra obtained specif- 
ically for constraining Ax/x with attached ThAr calibration lamps 
for each spectrum. This minimizes the systematic error of wave- 
length calibration. (3) Very high (~1 kms~'per channel) or high 
(2-4 kms~'per channel) resolution 21-cm spectra are used. (4) 
As we could get repeated observations for 21-cm absorptions we 
are able to identify the RFI related problems in the absorption pro- 
files. (5) We also estimate Ax/x by using only absorbers for which 
the background sources are unresolved even at milliarcsec scale in 
VLBA images. 

Although tight constraints on the variation of fundamental 
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Figure 14. Comparison of Aa/a estimated in this work with other mea- 
surements in literature. The dashed dotted line and the surrounded solid 
box show our measured Aa /a and its error and the filled square is the one 
from lSrianand et all 1201 Ol) based on the same method and assuming /.( and 
gp are not changing with time. 
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Figure 15. Comparison of A/i//j estimated in this work with other measure- 
ments in literature. The dashed dotted line and the surrounded solid box 
show our measured A/i/yii and its error based on our Ax/x measurements 
and assuming a and g p do not vary. Apart from the fi lled square at z~3.2 
which is calculated from Ax/x of Srianand et al. (2010) with the assumption 
of non-variation of other constants, the rest of the measurements for z > 2 
are based on the analysi s of H2 elec tronic transitions. The empty t riangl e 
towards up is from lKing et alj 1201 lh. empty star fromlMalec et all l2010h . 
empty square with a large bar from |y an Weerde nburg et alT i l2*oTlD . empty 
circle with a large bar from | Wendt & MolarolfcOl lb . filled triangles towards 
down from iThompson et alTl2*00^) . and filled triangles towards up from 
iReinhold et alj 120061) . The two measurements at z < 1 are based on the 
molecul ar inversion an d rotational tra nsitions. T he filled and empty squares 
are from Henkel et al. (2008) and Murphy et al. (2008) respectively. 



constants are obtained by comparing 21-cm and UV-optical red- 
shifts, the method is not exempt of systematics as all other methods 
in this field. As it is has already been discussed, the main source 
of uncertainty on Ax/x is related to the assumption used to asso- 
ciate one of the several UV-optical components to the 21-cm ab- 
sorption. By choosing absorbers with compact background radio 
sources at mas scale one can minimize the uncertainties related to 



the possibility that optical and radio sightlines are different. Dif- 
ferent methods have been imple mented to associate the UV-optical 
component to the 21-cm one. iTzanavaris et al.l d2007l) associated 
the pixel with strongest UV-optical optical depth to the pixel with 
strongest 21-cm optical depth. In this work we follow the same idea 
but using components of si multaneous Voigt profile fitting models. 
Using the same method as Tzanavaris et alj d2007l) . we find Ax/x 
= (3.6 ± 3.1) x 1CT 6 . Keeping this in mind we will now discuss the 
implication of our constraint on Ax/x on the variation of individual 
constants that constitute x. 

As x = g p a 2 /yu, its variation can be related to the variation of 
g p , a, and/or fi via Ax/x = Ag p /g v + 2 x Aa/a - Afi//i. Therefore 
the constancy of x can be related either to the constancy of all the 
three constants, g p , a and /i, or to some complicated combination 
of variations of these constants with an overall null effect on x. 

Assuming /j and g p are constants then our measured Ax/x 
translates to Aa/a = (0.0 ± 0.8) x 10~ 6 which is one of the most 
stringent constraint on the variation of a. In Fig[l4]we summarize 
the available constraints on Aa/a from the literature. It is clear that 
our measurements are consistent with Aa/a= (0.1 ± 1.5) x 10~ 6 
(and a factor 2 bet ter than that) found bv lSrianand et al.ld2007l) and 
with the results of IWebb et ail d201ll) f or z < 1.8 UVES data. E ven 
if we use the conservative approach of ITzanavaris et all d2007l) we 
get Aa/a = (1.8 ± 1.5) x 10~ 6 which is as good as the results from 
the MM method. 

IWebbet~aHj201ll) used a combined set of absorbers observed 
with VLT/UVES and Keck/HIRES to conjecture about the possible 
presence of a spatial dipole pattern in the variation of a. Their best 
fitted model indicates a spatial dipole in the direction with right 
ascension 17.5±0.9 h and declination -58 ± 9 deg, significant at 
the 4.2cr level with Aa/a = Arcos(0), where A = (1.1 ± 0.25) x 
10~ 6 GLyr _1 , r(z) = ct{z), c being the speed of light and t(z) the look 
back time, and is the angle on the sky between the quasar sight 
line and the best fit dipole position. We calculated the prediction of 
this m odel for our sample using the same cosmology as We bb et all 
J201 ]]). We also calculated the error on Aa/a using 



Aa 
<r( — ) 

a 



Aa/a\ , / Aa/a 



cos(0) 



(3) 



The measured and accordingly the predicted values for Aa/a 
are presented in columns 11 and 12 of Table [6] Last column of 
Table|6]shows the difference of Aa/a between the dipole and our 
measurement and its error. For the sight lines towards J 1623 +07 18 
and J0501-0159 the dipole predicts a variation of a at, respec- 
tively, the ~2.4<x and 2.9<r level. In both cases, our null results 
on the a-variation are inconsistent with this dipole prediction at 
more than 90% confidence level. For the other two sight lines to- 
wards J2340-0053 and J2358-1020 our measurements are con- 
sistent with null variati on as also predicted by the dipole. In 
the case of J1337+3152 ISrianand et all d201Cj|) measured Ax/x = 
-(1.7 + 1.7) x 10~ 6 for the absorber at z abs = 3.174 which translates 
to Aa/a = -(0.9 + 0.9) x 10~ 6 . The dipole prediction in this case is 
Aa/a = -(2.70 ± 1 .9) x 10~ s with a difference of -( 1 .80 ± 2. 1 ) x 10~ 6 
with the measurement. The difference between these five measure- 
ments and the dipole predictions results in a^ 2 of 14.5. The prob- 
ability of x 2 > 14.5 is ~ 1% which implies that the existence of 
a dipole is not favored by our measurements. More independent 
measurements especially towards systems where the dipole pre- 
dicts large variations will be useful to confirm/refute the existence 
of the a-dipole at higher significant level. 

Assuming that a and g p have been constant we derive Afi//j 
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= (0.0 ± 1.5) x 10~ 6 . Fig.fBlcompares our results with other direct 
measurements of Afi/fi obtained either using rotational transitions 
of H2 and HD molecules (for z > 2.0) or based on the comparison 
of NH 3 inversion transitions with some rotational transition lines 
(e.g. CO, CS, HC 3 N; for z < 1.0), While the constraints we get are 
not as good as the one obtained using NH 3 they are very stringent 
compared to those based on H 2 at z > 2. What is more interesting 
is that our measurements fill the redshift gap between NH 3 and Hi 
based measurements (see Fig.|15t. 

If we us e the lcr constraints on Aa/a found fo r z < 1.8 ab- 
sorbers (from lsrianand et al.l2007UWebb et al.l201 1) and Afi/fi es- 
timated at z ~ 0.7 using NH 3 dKanekaj|201ll) . considering they are 
valid at z ~ 1.3, we get Ag p /g p < 3.5 x 10~ 6 (l<x) from our Ax/x 
measurements. 

In summary, using 21 -cm and metal UV absorption lines we 
are able to derive stringent constraints on the variation of a, yu, 
and g p . As discussed before the best estimate on A/j/yU at z <1 is 
obtained by comparing the frequencies of NH 3 inversion transi- 
tions with rotational transitions of other molecules. The existing 
two measurements are towards the line of sight of two well known 
gravitationally lensed BL Lacs (B 0218+357 and PKS 1830-211) 
that show complex radio morphologies. As different transitions oc- 
cur at different frequencies the dependence of the background radio 
stru cture on frequency is an important source of systematic error 
(see iMurphv et ai]|2008l : lKanekarll201 lj). Therefore detecting NH 3 
and other molecules towards unlensed compact radio sources is im- 
portant to constrain A/i/fi. Unlike NH 3 and other complex heavy 
molecules, 21 -cm absorption is more frequently detected towards 
normal radio sources covering a wide redshift range. The main 
source of systematics in this method is related to how accurately 
the 21-cm absorption component is associated to the corresponding 
metal line component. More measurements towards compact radio 
sources are needed to address this issue adequately. Future blind 
searches for 21-cm absorption using the upcoming Square Kilome- 
ter Array (SKA) path finders hopefully will provide a large number 
of suitable targets to perform such measurements. 
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APPENDIX A: RESULTS OF CORRELATION ANALYSIS 
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Table Al. Shifts of individual spectra relative to the combined one in J0108-0037 at the position of different absorption profiles 



species 


EXP1 


EXP2 


EXP3 


EXP4 




Cms-') 


(ms- 1 ) 


(ms- 1 ) 


(ms- 1 ) 


(1) 


(2) 


(3) 


(4) 




Nin ,11454 


656± 505 


-826± 681 


531±458 


-536± 463 


Ni H/11467 


-131± 595 


-1081+ 835 


-973± 764 


484± 534 


Ni ii A 1502 


-1094± 556 


925± 887 


62± 816 


65± 941 


Ni nil 703 


-391+ 359 


169+ 473 


1306± 778 


161+ 618 


Ni H/11709 


-11± 184 


230±211 


-408± 227 


43± 195 


Ni H/11741 


-145± 142 


275± 164 


-160± 171 


-98± 162 


Ni H/11751 


-496± 204 


240± 221 


165+ 223 


-152± 211 


Cm ,12056 


-149± 109 


109± 121 


133+ 135 


-80± 121 


Cr n/12062 


-83± 117 


-168+ 130 


99± 139 


151± 132 


Crii/12066 


-66± 162 


15± 176 


-267± 202 


220± 179 


Mn ii ,12576 


32± 102 


115+ 111 


-44± 110 


-123+ 105 


Mnn ,12594 


-176+ 133 


-220± 151 


181± 150 


156+ 142 


Mn n,12606 


-187+ 149 


607± 173 


-181+ 168 


-140± 165 


Zn n,i2026 


-207± 166 


234± 188 


-209± 194 


170± 184 


Zn n,12062 


-108± 166 


79± 175 


-267± 198 


218+ 182 


Fen ,11611 


-282± 228 


74± 290 


752± 305 


-373± 267 


Feii/12249 


-145+ 103 


80± 110 


-118+ 116 


182+ 118 


Fe H/12260 


-358± 90 


300± 100 


-56± 105 


183+ 106 


Sin ,11808 


-237± 108 


132± 127 


239± 138 


-50± 124 


weighted mean 


-171 


+127 


-14 


+38 


weighted standard deviation 


152 


207 


223 


166 


weighted standard deviation of all exposures 




134 







Table A2. Shifts of individual spectra relative to the combined one in J 1623 +07 18 at the position of different absorption profiles 



species 


EXP1 


EXP2 


EXP3 


EXP4 




(ms- 1 ) 


(ms- 1 ) 


(ms- 1 ) 


(ms- 1 ) 


Ni ii ,11454 


2075 ± 573 


-414± 582 


94± 393 


-840± 364 


Ni ii A 1703 


2230±1013 


-531+ 682 


378± 446 


-414±419 


Ni iiil709 


-161+ 465 


520± 414 


485 ± 499 


-1040± 414 


Crii/12062 


177± 564 


-180+ 341 


-191+ 530 


242± 274 


Cm ,12066 


-434± 687 


461± 633 


335± 745 


-684± 602 


Znn ,12062 


-1434± 643 


461 ± 633 


72± 615 


460± 408 


Mn iii2576 


498± 426 


0± 346 


235 ± 420 


-217+ 224 


Mn n,12594 


-46± 458 


80± 303 


572± 386 


-185+ 213 


Mn iii2606 


475 ±481 


-1068+ 735 


-1073± 574 


194± 283 


Fen ,11611 


2002± 612 


1394+ 770 


766±1559 


-432± 335 


Fen ,12249 


-555± 389 


500± 273 


-415+ 424 


-70± 305 


Fe iii2260 


342± 320 


635± 322 


-632± 366 


-164± 194 


Fe n/12374 


-55± 184 


678± 161 


-913+ 193 


8± 110 


Fen i2382 


-275± 142 


460± 139 


-329± 160 


25± 89 


Feii/12586 


196± 160 


451+ 139 


-237± 158 


-307± 87 


Sin A 1808 


142± 261 


224± 220 


414± 253 


-381+ 180 


weighted mean 


+44 


+401 


-233 


-145 


weighted standard deviation 


543 


314 


480 


251 


weighted standard deviation of all exposures 




313 
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Table A3. Shifts of individual spectra relative to the combined one in J2340-0053 at the position of different absorption profiles 



species 


EXP1 


EXP2 


EXP3 


EXP4 


EXP5 


EXP6 


HIRES 




(m s -1 ) 


(m s -1 ) 


(m s _1 ) 


(ms~') 


(ms~') 


(ms~') 


(ms~') 


Nin /1 1709 


-645± 358 


212+ 428 


144± 336 


759± 455 


380± 412 


74± 345 


2530± 296 


Nin ,11741 


-81± 380 


5± 361 


447± 296 


668± 321 


97± 282 


-958± 317 


626± 232 


Nin /tl751 


-409± 513 


-122± 392 


588± 297 


-180± 373 


-406± 354 


142± 332 


65± 521 


Cm ,12056 


-169+ 759 


303± 868 


139+ 761 


-307± 922 


263± 811 


-41± 555 


-43± 196 


Cm ,12062 


-206± 749 


184+ 747 


217± 574 


49± 927 


-213+ 574 


-467± 606 


244± 268 


Cm ,12066 


728±1026 


-646±1201 


233± 988 


-561±1246 


74±1179 


-446±1042 




Zn ii,12026 


155+ 425 


130+ 486 


-88± 377 


-692± 389 


205 ± 290 


194+ 325 





Mn ii,12576 


-197± 351 


219+ 407 


487± 400 


-10± 384 


-13± 368 


-209± 290 




Mmi ,12594 


-62 1± 375 


129± 370 


122± 384 


-41 ± 363 


476± 361 


-110+ 274 




Mnn ,12606 


-143+ 317 


-439± 405 


885± 310 


-139+ 361 


436± 298 


-1199+ 284 




Fen ,11608 


-309± 66 


117+ 72 


275± 62 


-48± 81 


113+ 68 


-205 ± 65 


47±48 


Fen/11611 


-15± 408 


128+ 385 


66± 287 


281 ±493 


-325± 532 


-764± 456 


227 1± 539 


Fen ,12249 


-103± 184 


-104± 260 


-221+ 198 


-307± 195 


171+ 183 


138+ 149 


-433± 193 


Fen ,12260 


-1272± 138 


208± 133 


110+ 117 


43± 136 


594± 109 


132+ 101 


-27± 21 


Fen ,12374 


-308± 41 


185± 41 


230± 37 


35±40 


315+ 35 


-294± 32 


148± 57 


Sin A 1808 


-380± 139 


-137+ 160 


357± 144 


122± 154 


249± 157 


-165+ 134 


-229± 92 


weighted mean 


-352 


146 


236 


18 


283 


-238 


110 


weighted standard deviation 


296 


127 


148 


173 


184 


238 


438 


weighted standard deviation of all exposures 






204 


k 









★standard deviation in HIRES column is not included in the averaged standard deviation. 



Table A4. Shifts of individual spectra relative to the combined one in J2358-1020 at the position of different absorption profiles 



species 


EXP1 


EXP2 


EXP3 


EXP4 


EXP5 


EXP6 


EXP7 




(ms- 1 ) 


(ms" 1 ) 


(m s" 1 ) 


(ms" 1 ) 


(ms- 1 ) 


(ms" 1 ) 


(ms" 1 ) 


Zn ii ,12026 


-509± 549 


-477± 376 


-744± 629 


1083± 395 


173+ 266 


165+ 326 


-511+ 553 


Mn ii ,12576 


-555± 521 


73±419 


-658± 420 


-488± 406 


150+ 317 


322±413 


1316+ 578 


Fe ii ,12249 


-13±747 


-317+1013 


-370± 596 


-203± 452 


753± 730 


397± 456 


-568± 473 


Fe ii ,12344 


-19± 234 


88±251 


-257± 247 


-311+ 189 


32± 151 


29()± 201 


72± 201 


Fen ,12374 


-800± 294 


- 118+ 213 


-35±261 


10± 177 


205± 156 


188+ 174 


392± 264 


Fen ,12586 


-130+ 268 


331+ 243 


24± 245 


-213+ 166 


139+ 165 


258± 180 


-178+ 246 


Si nil 808 


-20± 577 


-619+ 661 


-1897± 719 


-497± 499 


241±421 


612± 527 


1975± 768 


weighted mean 


-275 


+70 


-242 


-129 


+142 


+ 147 


+ 111 


weighted standard deviation 


357 


288 


431 


374 


113 


260 


561 



weighted standard deviation of all exposures 170 



